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A robust control law design methodology is presented to stabilize the X-56A model and 
command its wing shape. The X-56A was purposely designed to experience flutter modes in 
its flight envelope. The methodology introduces three phases: the controller design phase, 
the modal filter design phase, and the reference signal design phase. A mu-optimal controller 
is designed and made robust to speed and parameter variations. A conversion technique is 
presented for generating sensor strain modes from sensor deformation mode shapes. The 
sensor modes are utilized for modal filtering and simulating fiber optic sensors for feedback 
to the controller. To generate appropriate virtual deformation reference signals, rigid-body 
corrections are introduced to the deformation mode shapes. After successful completion of 
the phases, virtual deformation control is demonstrated. The wing is deformed and it is 
shown that angle-of-attack changes occur which could potentially be used to an advantage. 
The X-56A program must demonstrate active flutter suppression. It is shown that the virtual 
deformation controller can achieve active flutter suppression on the X-56A simulation 
model. 
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D a = frequency weighting matrix 

d re f = reference deformations 

d(x C) y c , z c , t) = vector of deformations at positions x c , y c , z c at time t 

E = rational approximation matrix of lag terms 

FEM = finite element model 

FOS = fiber optic sensors 

F t = objective function for the i th mode 

/ = number of aerodynamic lag states 

freq = frequency 

G = plant 

GAF = generalized aerodynamic force 

G v = set of plants over a velocity range 

Gdes = design plant 

g = damping 

I = identity matrix 

i A = input uncertainty signal 

i = imaginary number 

J = generalized damping matrix 

K = controller 

k m = stability margin 

L = characteristic length 

LFT = linear fractional transformation 

LWLE = left-wing leading edge 

LWTE = left-wing trailing edge 

l r = sensor row index vector for reference deformations 

M = generalized mass matrix 

M a = matrix transfer function 

m = number of modes in model 

m r = modal column index vector for reference modes 

N = number of nodes in finite element model 

NS = nominal stability 

NP = nominal performance 

N a = matrix transfer function 

n = number of states 

7i r = rigid body sensor noise 

n s = fiber optic sensor noise 

n u = input control noise 

0 S = original sensor position before modal deformation 

Os = original sensor position after modal deformation 

OP = line segment distance from point 0 to point P 

O'P' = line segment distance from point O' to point P' 

o = number of measurements 

o A = output uncertainty signal 

P = sensor location point 

P g = generalized plant 

P + = sensor position on the fore side of original sensor before modal deformation 

P+ = sensor position on the fore side of original sensor after modal deformation 

P_ = sensor position on the aft side of original sensor before modal deformation 

Pi = sensor position on the aft side of original sensor after modal deformation 

P ro Ju(o^pl)OsP+ - projection of directed segment on unit vector in the direction of 0 S P + 

P ro Ju(o^pl)OsPl - projection of directed segment 0' S P1 on unit vector in the direction of 0 S P_ 

P r °ju(dP)0'P' = P r °j ec ti° n °f directed segment O'P' on unit vector in the direction of OP 

Q = generalized aerodynamic force matrix 
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Q(s) = rational function approximation of the generalized aerodynamic forces 

q(t) = vector of modal displacements at time t 

q (t) = vector of modal velocities at time t 

q (t) = vector of modal accelerations at time t 

£/;(0 = i th modal displacement defined at time t 

q re f = modal reference displacement command 

R = diagonal matrix of lag roots 

RP = robust performance 

RS = robust stability 

RWLE = right-wing leading edge 

RWTE = right-wing leading edge 

(Ore/ = reference command of the argument 

5 = generalized stiffness matrix 

SFOS = simulated fiber optic sensors 

SP = short-period mode 

SWBT = symmetric wing bending and torsion flutter mode 

SW1B = symmetric wing first bending 

SW1T = symmetric wing first torsion 

Sj = input sensitivity matrix 

S 0 = output sensitivity matrix 

S d - strain-displacement transformation matrix 

5 = Laplace variable 

sim = simulation 

s m ( x c , y c ,z c ,t ) = strain at measurement locations at time t 

7(0 = right-hand-rotation matrix 

Tj = input complementary sensitivity matrix 

T 0 = output complementary sensitivity matrix 

t = time 

u = control signal vector 

u(t) = control input signal vector at time t 

u(-) = unit vector 

V-f = frequency versus velocity 

V-g = damping versus velocity 

Vdes = design speed 

Voo = free- stream velocity 

v = error vector from reference input minus measurements 

W = weighting 

WF = wing flap 

WF1R = wing flap first right 

WF1L = wing flap first left 

WF2R = wing flap second right 

WF2L = wing flap second left 

WF3R = wing flap third right 

WF3L = wing flap third left 

WF4R = wing flap fourth right 

WF4L = wing flap fourth left 

W A = matrix of proper additive weights 

W I = matrix of proper input weights 

W 0 = matrix of proper output weights 

W s = matrix of proper sensitivity weights 

W T = matrix of proper complementary sensitivity weights 

W u = matrix of proper control output weights 

w = disturbance inputs 

X = data matrix 

x(t) = state vector at time t 
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state derivative vector at time t 
Cartesian coordinate in the x-direction 
vector of airframe states 
reference airframe states 
vector of airframe states at time t 

vector of actuator displacements, velocities and accelerations at time t 
c.g. location of the aircraft in the x c direction 
vector of modal displacements 

uni-axial coordinate in a moving frame always in the direction of the SFOS layout 

vector of aero lag states at time t 

initial state at time 0 

modal displacements 

vector of modal displacements at time t 

vector of modal velocities at time t 

un-deformed location of the aircraft nodes in the x c direction 

deformed location of aircraft nodes for the i th mode shape in x c 

estimated modal displacement states 

rational approximation matrix of lag terms 

output sensor vector 

output sensor vector at time t 

Cartesian coordinate in the y-direction 

c.g. location of the aircraft in the y c direction 

un-deformed location of the aircraft nodes in the y c direction 

deformed location of aircraft nodes for the i th mode shape in y c 

number of actuator states 

Cartesian coordinate in the z-direction 

goal regulated KS states 

goal regulated sensitivity states 

goal regulated complementary sensitivity states 

c.g. location of the aircraft in the z c direction 

un-deformed location of the aircraft nodes in the z c direction 

deformed location of aircraft nodes for the i th mode shape in z c 

angle of attack 

angle of sideslip 

structured uncertainty 

input uncertainty 

change of velocity 

change of angle of attack 

change in pitch angle 

scalar global displacement in x c 

forward change of aircraft’s center of mass from trim 

scalar global displacement in y c 

side change of aircraft’s center of mass from trim 

scalar global displacement in z c 

heave change of the aircraft’s center of mass from trim 

scalar global angle about y c axis 

scalar global angle about x c axis 

scalar global angle about z c axis 

change in deformation at LWLE 

change in deformation at LWTE 

change in deformation at RWLE 

change in deformation at RWTE 

change in SW1T modal displacement 

change in SW1B modal displacement 
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axial strain measured on aft side of a sensor 
pitch angle 

structured singular value 
free- stream density 
maximum singular value 

deformation modal matrix, a collection of mode shapes, <p m 
bank angle 

i ih natural deformation mode shape 

= i th natural deformation mode shape defined over x c ,y c , z c 
strain mode matrix, a collection of strain mode shapes 
strain matrix defined at the SFOS locations 
yaw angle 
i th strain mode 
m th strain mode 
frequency 

derivative in axial beam deformation in coordinate direction x x 
Hardy space norm 


I. Introduction 

T he NASA Aeronautics Research Mission Directorate (ARMD) is focusing on developing technology for 
dramatically reducing noise and emissions while improving performance for transport aircraft. 1 New technology 
must support the decrease of aircraft fuel and energy consumption by 60 percent within a timeline of 15 years. One 
method of achieving this goal is to design new lightweight aircraft. The consequence of using a lighter structure with 
increased flexibility is that the structure is more likely to experience adverse aeroelastic effects such as flutter and 
divergence. 2 Researchers suggest that flexible aircraft may be more susceptible than traditional aircraft to gusts and 
maneuver loading. 3 Aircraft designed with lightweight flexible structures will require active control technology to be 
feasible alternatives to stiffer aircraft. 4 

Research into active flutter suppression (AFS) technology on flexible aircraft is ongoing. The United States Air 
Force Research Laboratory (AFRL) Multi-utility Aeroelastic Demonstration (MAD) Program has developed the 
Multi Utility Technology Test-bed (MUTT), also known as the X-56A (Lockheed Martin, Bethesda, Maryland). 5 * 
The X-56A is a flutter demonstrator designed for the sole purpose of testing various AFS solutions. Lockheed 
Martin is the prime contractor for design and development of the vehicle. The X-56A stems from previous work 
between Lockheed Martin and the AFRL to design and develop high-altitude, subsonic, long -endurance autonomous 
aircraft. 6,7 

The X-56A is a remotely -piloted aircraft with a stiff center body and flexible detachable wings. Due to the high 
flexibility of the wings, at certain speeds within the flight envelope the aircraft modes couple unfavorably, leading to 
three flutter instabilities. A control system must be developed with existing sensors and actuators to suppress the 
instabilities. Robust control of these flutter modes will demonstrate that lightweight structures can be implemented 
and flown safely. As part of its sensor array, the X-56A aircraft wings will be instrumented with fiber optic sensors 
(FOS) and fiber Bragg gratings. 8 Fiber optic sensors fall into a class of distributed sensors and can measure strain at 
thousands of locations. A research objective in the X-56A program is to use these sensors in a feedback control 
system. This is the problem this paper specifically addresses. 

Active flutter suppression technology may be extended to not only suppress flutter modes, but to also control the 
shape of the aircraft. Research has been published which indicates that drag can be reduced due to structural 
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deformation changes. 9 The X-53 Active Aeroelastic Wing aircraft (The Boeing Company, Chicago, Illinois) warped 
its wing to achieve higher roll rates, improving its maneuverability. 10 Research in shape control is ongoing. Smart 
technologies such as adaptive control surfaces, 11 adaptive ribs, 12 piezoelectric, 13 and piezoceramic actuators 14 are 
under investigation. Smart actuation will not be available on the X-56A, so this paper investigates the use of 
trailing-edge control surfaces for shape control. 

Presented in this paper is a robust control law design methodology for stabilizing X-56A aircraft models and 
controlling the wing shapes. The models contain all six rigid-body degrees of freedom (DOF), flexible modes, and 
ten control surfaces and flight actuators. To use the FOS, the methodology requires implementation of a 
least-squares estimation (LSE) modal filter. In previous work, 15 the LSE modal filter was tested on a clamped 
aeroelastic plate model, reinforced with spars and ribs and controlled by two trailing-edge control surfaces. In that 
work, simulated shape data at the sensor locations was used for sensor inputs to the shape controller. Here, we 
utilize simulated fiber optic sensor (SFOS) measurements for the shape controller. We derive strain modes from the 
deformation mode shapes and use these modes for both modal estimation and strain simulation. 

The use of spatial filtering (that is, modal filtering) and distributed sensing enables practical shape control. Point 
sensors such as strain gages or accelerometers could be utilized to support shape control. But arrays of point sensors 
must be frequency-filtered or placed optimally to capture important modal coordinates of interest, or both. 16 
Distributed sensors with modal filtering do not introduce unwanted dynamics into the problem, and the modal 
displacement estimates will be resistant to the effects of residual modes . 17,15,18 

Three objectives are achieved in this paper. The first objective is to verify that the modal filter is a valid 
estimator for an aircraft. The second objective is to show that virtual deformation control (see Ref. 15) on an aircraft 
in flight is achievable. This is accomplished with SFOS feedback 8 and trailing-edge control surface effectors. The 
third objective is to show that a controller designed to incorporate a modal filter can be made to be robust to 
structured uncertainties and speed variations. We demonstrate that the popular /^-optimal approach 19,20 can be used 
to design a modal controller which achieves multiple-input multiple-output robust stability and performance. 


II. The X-56A Model 

The wings of the X-56A model are fitted with four trailing-edge control surfaces. There are two body flaps on 
the rigid center body. All of these control surfaces are available to the controller for both AFS and flight control. 
The medium-sized unmanned aerial vehicle aircraft will be operated in experimental studies at subsonic speeds and 
at low altitude. A top-down view of the X-56A aircraft is shown in Fig. 1. 



Figure 1. Top-down configuration of the X-56A (L: left; R: right) (Figure courtesy Lockheed Martin, 
Bethesda, Maryland). 
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To simulate fuel in the wings for a passenger transport aircraft, water tanks are included on each wing. This 
generally makes control design more difficult because the controller will either need to adapt or be robust to the 
varying modal properties. This challenge may be compared to the wing stores AFS problem of the Northrop YF-17 
airplane. 21 For simplicity, the simulation model used here assumes that the water tanks are empty. The aircraft has 
no rudder, so its yaw control will be limited. To yaw, the control surfaces must move in opposing directions to 
generate more drag on one side of the vehicle. To identify the flutter characteristics of the vehicle, a normal modes 
analysis is first completed. 

A. Normal Modes Analysis 

The normal modes analysis was completed in MSC Nastran 22 (MSC Software Corporation, Santa Ana, 
California) on a detailed finite element model (FEM) developed by Lockheed Martin. The FEM corresponding to 
the water-empty case was utilized. The detailed FEM is shown in Fig. 2. 



Figure 1. The finite element model of the X-56A aircraft. 

The FEM contains over 6,000 nodes. The FEM has not yet been updated through a ground vibration test. For model 
development and flutter analysis, the modes were computed in MSC Nastran using the Lanczos eigensolver. The 
modes are shown here amplified by a large amount (for visualization) in Fig 3. 
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Figure 3. Normal mode shapes: a) SW1B; b) AW1B; c) SW1T; and d) A WIT. 
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Only four modes are shown here, but 14 flexible modes were included in all of the X-56A models. Fourteen 
were found to be sufficient to capture the salient flexible motion characteristics. Visual inspection reveals that the 
first two mode shapes appear to be strongly coupled with rigid-body motion. The coupling is pronounced in the 
symmetric wing first bending (SW1B) mode, where rigid-body pitch and heave are observed. The antisymmetric 
wing first bending (AW IB) mode has a substantial rigid-body roll component. The symmetric first wing torsion 
(SW1T) mode has a slight pitch coupling which is difficult to see. The antisymmetric wing first torsion (AW1T) is 
slightly coupled with roll. The flutter analysis presented below elucidates how the modes couple in flight. 

B. Flutter Analysis 

To determine the theoretical flutter margins of a vehicle the eigensolution of the aeroelastic system equation is 
solved 2 as shown in Eq. (1): 


Mq(t) +]q(t) + Sq(t ) = -p m V^Qq(t) 


( 1 ) 


where M E M mxrn is the generalized mass matrix, / E M mxrn is the generalized damping matrix, S E M mxrn is the 
generalized stiffness matrix, Q E IR mxm is a generalized aerodynamic force (GAF) matrix, and q(t) E M mxl are the 
modal displacements of the vehicle. The aerodynamic interaction of the modes with the structure is captured by the 
GAF matrices. For flutter analysis, the GAFs must be converted from individual frequencies to Laplace domain. To 
achieve this, the GAF matrices, Q , are computed over a range of reduced frequencies and are fit in a rational 
function approximation (RFA) as shown in Eq.(2): 


L L 2 

Q(s) = A a + ttA-lS + — 2 A 2 s 2 + Y 

too too 


L 


Es 


( 2 ) 


where A 0 , A ± , and A 2 are rational approximation matrices, Y and E are rational approximation matrices of lag terms, 
R is a diagonal matrix of lag roots, L is the characteristic length, and is the freestream velocity. Roger’s method 23 
with three aerodynamic lag orders was used for the aerodynamic approximation in the ZONA Technology Inc. 
(Scottsdale, Arizona) ZAERO software package. 24 After substituting Eq. (2) into Eq. (1), and performing an inverse 
Laplace transform, the time domain aerodynamics in the equations of motion is achieved. It is important to verify 
that the frequency domain aerodynamics match the time domain aerodynamics for model accuracy. Therefore, in the 
flutter analysis, both frequency and time domain solutions are presented together in Fig. 4. 
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SP(freq) — 0 — SWIB(freq) — * — AWIB(freq) — * — SWIT(freq) — * — AWIT(freq) 


H SP (time) — O — SWIB(time) — * — AWIB(time) — ☆ — SWIT(time) — x — AWIT(time) 



Figure 4. Flutter analysis: a) velocity versus frequency, V-f ; and b) velocity versus damping, V-g. 
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At first glance, Fig. 4 indicates that the aerodynamic approximation for the state space models matches the 
frequency domain aerodynamic solutions. The primary purpose of Fig. 4 is to reveal the flutter characteristics of the 
X-56A model. Theoretically, flutter appears at a velocity in the model when the damping, ( g ), goes to 0. In practice, 
flutter predicted in this way is conservative, as structural damping is unaccounted for. 

Figure 4(b) reveals that three flutter modes exist. The first mode that goes unstable (that is, crosses 0 g) is the 
body freedom flutter mode (BFF) at a ratio of design speed, V des , of approximately 0.82. The BFF frequency is 2.7 
Hz (16.9 rad/s). The viscous damping ratio for the BFF is -0.13. Recall that the g value is equal to the negative of 
twice the viscous damping. 25 Therefore, positive g corresponds to an unstable flight condition. The symmetric wing 
bending and torsion flutter mode (SWBT) appears at approximately 89 percent of design speed at 6.5Hz (41 rad/s). 
The viscous damping ratio of the SWBT is lower than the BFF at -0.06. The velocity versus frequency plot, V-f, 
indicates that unfavorable coupling occurs between the SW1B and SW1T modes, similar to what one might expect 
when performing a flutter analysis on a clamped plate. 25 Typically, the coupled modes’ frequencies shift down but 
do not coalesce. The antisymmetric wing bending and torsion flutter mode (AWBT) appears at 94 percent of design 
speed at 4.9 Hz (30.7 rad/s). The damping ratio of the AWBT is -0.024. This represents another unfavorable 
coupling between the AW IB and AW IT modes. 

Figure 4 shows three strong flutter modes in the flight envelope which must be designed for in the active control 
system. The viscous damping of each mode decreases (becoming more unstable) very quickly, relative to velocity 
changes. This is especially true for the BFF. The flutter characteristics (damping ratios and frequencies) presented 
above represent the salient characteristics of the state space models over a range of speed. The development of the 
X-56A state space models is further discussed below. 

C. X-56A State Space Modeling 

The X-56A aeroservoelastic state space model assumes the matrix time derivative form for simulation and 
control design shown in Eq. (3) (see Ref. 24 or 26): 

x(t) = Ax(t) + Bu(t ) 

y(t) = Cx(t ) + Du(t ) W 0 ^ ’ 


with the initial state x(0) = x 0 . The n- dimensional vector x(t) is referred to as a state vector and at any discrete 
time during a simulation can be accessed to give the current “state” of the system. The o-dimensional vector y is the 
system measurements. According to Eq.(3), only the current state and the a-dimensional input u is required to know 
the state in the next time step. The A, B, C and D matrices are real constant matrices with nxn, nxa, oxn and 
o X a dimensions. States are not unique to any given system; however, the X-56A state vector is defined as shown 
in Eq. (4): 


*00 = {XAp(f). X e (t), x e (t), X lag (t), X act (t)} T 


( 4 ) 


where x AF (t ) E M^ xl is a vector of airframe states, x e (i) E R mxl is a vector of modal displacements, x^(t) E 
R mxl is a vector of modal velocities, Xi ag (t) E IR^ xl is a vector of aerodynamic lag states and x act (t) E R zxl is 
formed from a vector of actuator accelerations, velocities, and displacements. For the X-56A model, the state space 
models were calculated at small velocity increments over a range of speed. The resulting models are collected into 
the aggregate model, G v , to be used later for analysis. 

Figure 4 predicts that the design plant G des will be subject to three flutter modes. To verify that these modes are 
present at the design speed, the X-56A state space model was perturbed with a unit-scaled control deflection 
command to the right-wing control surfaces (see Fig. 1), lasting 0.001 s. The pertinent state space variations are 
plotted in Fig. 5. 
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The BFF in Fig. 5(a) is demonstrated by the interaction of the scaled states, rigid-body heave state, 8z cm , and modal 
displacement 8SWB1. The effect of gravity is also modeled. Note that altitude is lost over time. It is observed that 
the amplitude of the interacting states increases with time and oscillates at the same frequency, out of phase. The 
SWBT also illustrates the unfavorable coupling of modal displacements: 8SW1B and 8SW1T (see Fig. 5[b]). The 
nature of the coupling is difficult to discern because it appears random in nature. The AWBT shows an in-phase 
interaction of modal displacements: 8 AW IB and 8 AW IT. The scaled modal amplitudes are small but clearly grow 
in time, verifying that this flutter mode also exists at the design flight condition. Without control, the aircraft is 
clearly going unstable and will require AFS at a minimum. 

III. Methodology 

The model of the unstable X-56A aircraft was introduced above. Three flutter modes were identified, including 
BFF, SWBT and AWBT before the design speed. Active flutter suppression is desired to suppress the flutter modes; 
however, a higher objective of this work is to sense and command the shape of the aircraft. The available sensors 
include SFOS and flight sensors. It is shown how the SFOS may be used to estimate the modal displacements with 
the modal filter. This is followed by incorporating the modal filter into an aircraft shape control design paradigm. 
The overall methodology is designed to support the simulation framework given in Fig. 6. 


n u 



Figure 6. The virtual deformation control architecture for the X-56A model. 


In this simulation framework, there are three important steps which must be taken. First, the controller, K , must be 
designed. Second, the modal filter and SFOS must be developed; and third, the modal transformation must also be 
developed, which is used to generate modal displacement reference signals for the shape controller. The 
methodology below is presented, which supports the simulation shown in Fig. 6. The connections and variables are 
more completely described farther on in this paper. 

A. Modal Filtering in the Aircraft State Space Model 

The purpose of using the modal filter (see Fig. 6) is to estimate the modal displacements of a flexible structure. 
The modal filter was originally developed for the purpose of developing independent modal control systems, and 
was reliably found to remove the possibility of exciting residual modes in the control system. 17 In aircraft, 
frequency-based filters have been traditionally used for this purpose. 4 The first modal filter was derived using the 
modal expansion theorem, 27 which states that for structures with distinct natural frequencies, the deformation of all 
nodes of the structure d(x c ,y c ,z c ,t ) can be set equal to a linear combination of the natural mode shapes of the 
structure {<Pi(x c ,y c ,z c ), i = 1 ... N} as shown in Eq. (5): 


d(x c ,y c ,z c , t) 


^ qi(t)<pi(x c ,y c ,z c ) 


( 5 ) 
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The mode shapes are made orthogonal to the mass matrix, and are thus orthogonal to each other. The orthogonality 
property makes residual modes invisible to the estimator. Thus, it achieves spatially-filtered estimates without 
corruption from residual modes if enough sensors are available. 18 The highly-distributed nature of the SFOS is why 
the SFOS so naturally support modal filtering. Strain has a similar relationship, which has been used to define the 
strain-based modal filter. The relationship of the strain-based modal filter to the modal displacement state vector 
x e (t) is defined in Eq. (6) (see Refs. 28 and 15): 

x e (t) = ( W Tl V)- ll ¥ T s m (x Cl y Cl z Ct t ) + £ (6) 

where 4* = {S d 0;(x c ,y c ,z c ), i = 1 ... N] is a N X N strain matrix with N strain modes, s m (x c ,y c ,z c ,t) E R Nxl is 
the strain and 8 is the error. This represents the modal filter in its pseudo-inverse form. 28 Since the modal filter gives 
a partial state estimate of the full state vector, x(t), the form of the output matrix is simply identity for measured 
modes. 

The sensor output matrix, C, is a matrix of row vectors relating the sensory information to the state vector, x(t). 
Previously (see Ref. 15) the output matrix was developed for the case of a rectangular wing, where rigid -body 
motion was restrained. Here, the output matrix is adjusted for rigid-body state feedback concurrent with modal 
deformation state feedback. It is assumed that only m modes are retained for measurement. Similar to before, the 
output matrix is formed assuming all modal displacements, x e , and rigid-body states, x AF , are measured, as shown 
in Eq.(7). 


C A 


’ l ASX AS 
®mxAS 


®ASxm ®ASxm ®ASxf ®ASxz 
Imxm ®mxm ®mxf ®mxz. 


(7) 


Measuring all flexible states (or rigid states) is not required to adequately sense the vehicle state. Modes which are 
unstable and contribute strongly to the deformation of the structure are good candidates. Higher-order modes 
generally do not significantly contribute strongly to the deformation. 16,17 These modes may be cautiously ignored to 
reduce the size of the C matrix. An excellent method for choosing which modal displacements to feed back is 
identifying interacting modes from the flutter analysis presented in Fig. 4. Modes which interact strongly and lead to 
flutter are important to the feedback control system. 15 One may also look at the modal mass of each mode to 
determine which modal displacements will be significant to the response. 15 A convergence study 29 may be useful for 
choosing which modes to feed back to control a fluttering structure. 

B. Incorporation of Modal Filtering into Aircraft Control System Design 

Described above was how the modal filter may be incorporated mathematically into the state space model, which 
is in turn used for control design. The modal filter implementation was also demonstrated in Fig. 6. The 
implementation of the modal filter is split into phases. The phases suggested here include: the control design phase, 
the modal filter design phase, and the shape reference signal design phase. A methodology for incorporating a modal 
filter was given in Ref. 17, however, it was not specific for an aircraft or for shape control. The phases are presented 
as follows: 

7. Control Design: Phase I 

1) Define the robustness and performance requirements. 

2) Identify which modes are significant to the response of the model, through the flutter analysis (see Fig. 4). 

3) Update the output matrix to measure the identified modes in the state space model as in Eq. (7). 

4) Scale and reduce the order of the plant. 

5) Design the controller, , with the chosen control synthesis technique. 

6) Reduce the controller order. 

7) Verify that the performance and robustness requirements were met, while simulating the controller on a 
range of plants around the design plant. 

8) Iterate steps 2 through 7. 
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2. Modal Filter Design: Phase II 

1) After selecting the modes for feedback, place the sensors either using an optimal sensor placement 30 
technique or intuition of which modes are being measured. 

2) Form the sensor strain modal matrix, . 

3) Select an operator to calculate the modal displacements with the modal matrix. 

3. Shape Reference Signal Design: Phase III 

1) Identify the reference points on the vehicle for shape tracking, using the optimal sensor placement 
techniques or intuition. 

2) Correct the elastic mode shapes to remove the rigid-body translations or rotations, if any. 

3) Form a modal set-point transformation matrix with the corrected mode shapes. 

4) Verify that the flexible modes included in the controller lead to acceptable deformation control. If they do 
not, include more modes in the controller for integral tracking within the actuator bandwidth and repeat the 
control design phase. 

Depending on the platform in which the modal filter is implemented, some of the steps may not need to be taken. 
For example, if a plant is small enough, order reduction may not be required. During reference signal construction 
for a less-flexible structure, flexible modes may not be coupled with rigid-body translation and rotation. These 
modes will not need correction. 


IV. Control Design 

The methodology defines the way the modal filter fits into an aircraft shape control design paradigm. It was 
shown that only simple changes were required for the state space matrices to utilize the modal filter. See Eqs. (6) 
and (7). To design the controller as stated in Phase 1, the actual fiber optic sensor locations need not be determined 
yet; only the modal displacements which are desired to be controlled should be known. The requirements of the 
vehicle controller are now derived, and a p -optimal controller is designed, which can track both rigid-body states 
and the first two modal displacements. 

A. Control Design Overview 

It is desired to suppress the flutter modes (see Figs. 4 and 5) and track the rigid-body commands of the X-56A 
aircraft. In addition, one of the primary objectives is to demonstrate virtual deformation control on a flight vehicle. 
To this end, the contributing structural modal displacements will be tracked. In this regard, the shape of the aircraft 
will be controlled. The following discussion will be primarily concerned with how Phase I was accomplished. The 
first step of Phase I is to define the robustness and performance requirements. 

Since flutter is a potentially destructive phenomenon, the modal controller must be robust to uncertainty. This is 
especially important for a modal controller, since during a ground vibration test, mode shapes are typically not 
predicted exactly. 16 In fact, cross-coupling is a primary argument against using a modal filter for control. 
Cross-coupling can occur when any one mode shape in the modal matrix is not accurate and propagates into other 
modal displacement estimates through projection. 31 To reduce these uncertainties, the following uncertainty 
requirements were defined: 10 percent multiplicative uncertainty on the inputs and outputs, and 10 percent additive 
uncertainty on the scaled plant. Since flutter speeds are hard to predict precisely and structural damping is not 
modeled, the controller requirements must also be satisfied at off-design conditions. A notional requirement is that 
the controller must be robust to a 3 -percent speed variation. Since the flutter mode damping decreases so rapidly 
with speed variations (see Fig. 4), this requirement is seen as practical. 

Robustness to parameter uncertainty is highly desirable, but the controller must also meet performance 
specifications. The tracked measurements must respond to doublet inputs with low rise times and small overshoots. 
The performance must also be achievable in the aforementioned uncertainty conditions. The uncertain plant and 
required control system is summarized in Fig. 7. 
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Figure 7. The uncertain plant and the required controller. 

With the uncertainties in Fig. 7 in mind an ^-optimal control approach was taken. 19 The design approach uses a 
hybrid of performance weights and uncertainty weights. Robust stability (RS) is achieved if and only if a system is 
stable for all perturbed plants about the nominal model up to the worst-case uncertainty. 19 The robust stability 
condition which must be met for the ^-optimal controller is shown in Eq. (8). 


RS <=> i u(M A (im)) < l,Va) 


( 8 ) 


The ii is calculated over the frequency range with the relation shown in Eq. (9) (see Ref. 19): 

= min{k m \det(l - k M taO« = 0 for structured A. W) < 1) 00 

where k m is the stability margin, d is the maximum singular value, and A is the structured uncertainty. The transfer 
function matrix from the input of the uncertainty blocks to the outputs of them as shown in Fig. 7 is presented in Eq. 
( 10 ): 


-W,T, 


M a = 


W A S, 
,W 0 GS ! 


—W,KS 0 

-W a KS 0 

W 0 S 0 


—W,KS 0 

-w a ks 0 

~W 0 T 0 


, GEGy 


( 10 ) 


where Wj is a matrix of proper input weights, W A is a matrix of proper additive weights, W 0 is a matrix of proper 
output weights, and K is the controller. From this matrix, the salient sensitivity and complementary sensitivities 
from the M A structure are identified which correspond to the requirements. The magnitude of the singular values 
predict the performance of the control system. These closed-loop transfer functions are defined in Eq. (11) (see Ref. 
19): 


Tj = KGQ + KGY 1 
KS 0 = KQ + GKY 1 

S 0 = ( 1 + GKY 1 (11) 

Sj = (/ + KGY 1 
T 0 = GKQ + GKY 1 

where Tj is the input complementary sensitivity, S 0 is the output sensitivity, S 7 is the input sensitivity, and T 0 is the 
output complementary sensivitiy. To improve rejection of control input uncertainty, the controller was designed 
around the plant as shown in Fig. 8. 
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Figure 8. The control design framework. 

The design plant, G des , at the design speed, V des , was of order 130. To improve controller synthesis, G des was 
scaled by the full range of actuator movement and expected sensor output changes (see Ref. 19). The translational 
states, 8x cm , 8y cm , and 8z cm , and the velocity state were removed from the model. With the resulting 126th-order 
model, balanced reduction 32 was performed to bring the model order to 90 states. The selected states to be tracked 
were: SW1B and SW1T modal displacements, pitch angle, 6 , and bank angle, </>. The angle of attack a , angle of 
sideslip /?, yaw angle, 0 , and a body accelerometer were also sensed but were chosen to be suppressed. 

Traditional proper weights from a mixed K synthesis were utilized, along with multiplicative uncertainty at the 
plant inputs. The input uncertainty weight, Wj , was adjusted to achieve maximum amplitude near the actuator break 
frequency. Sensitivity weight, W s , was adjusted for integral tracking on tracking states and for suppression on 
suppression states. It was found that the break frequencies of the modal displacement performance weights had to be 
increased 10 rad/s relative to the airframe state weight break frequencies of 1 rad/s. The break frequency of the 
control weight, W u , was adjusted to 5 rad/s for reduced control surface movement. The break frequency of the 
complementary sensitivity weight, W T , was set to 30 rad/s to improve high frequency noise rejection. The 
uncertainty transfer function, N A , was calculated (see Fig. 8) to be as shown in Eq. (12): 


' W I T l -W^Sq 
-W sGSj -W s S 0 
~WJ, -W U KS 0 
_ W t GSj -W t T 0 _ 


( 12 ) 


One may verify that all of the pertinent closed-loop transfer functions corresponding to those shown in Eq. (11) are 
present in Eq. (12). Thus, by reducing the norm of N A , the robustness will increase and it will be easier to meet 
the specified requirements (see Fig. 7). For this control design architecture, the jU-optimal controller was computed 
using MATLAB’s Robust Control Tool Box. To find the controller, DK-iteration (See Ref. 19) was performed, 
which solves the iterative optimization problem shown in Eq. (13): 

mm(min||D ft ,JV 4 (J0i)«- 1 |L) ( 13 > 

The DK-iteration resulted in a 162th-order controller after some trial and error with the weights in Fig. 8. The 
controller was then internally balanced 32 and truncated to 44 states without a substantial loss of robustness or 
performance. This resulted in an 0~C o 0 norm of 3.29. It was difficult to meet both performance and robustness 
requirements with the scaled plant. Therefore, the desired Tfoo norm of 1 was not achieved. Rescaling may improve 
the controller. 
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Nominal stability (NS) was verified for the reduced-order controller by verifying that eigenvalues of the 
closed-loop system had real negative parts. An analysis of the singular values of the open-loop plant and controller 
as well as the closed-loop sensitivity functions given in Eq. (11) are presented in Fig. 9. 



Figure 9. Maximum singular values of open-loop (GK) and closed-loop sensitivity functions. 

At first glance, Fig. 9 reveals that most of the closed-loop transfer functions stay at or beneath 10 dB, 
corresponding to a magnification of the inputs to these transfer functions of approximately 3.2. The controller 
bandwidth is 25 rad/s, which is below the actuator break frequency. After this frequency, the complementary 
sensitivities fall sharply at 60 dB/decade. Therefore, noise will be attenuated well past 25 rad/s. 

Recall that it is desired that the controller be robust to uncertainty and, of course, disturbances. Robustness to 
input multiplicative uncertainty is strongly impacted by the maximum singular values of Tj (see Eq. [10]). The 
maximum singular value of Tj within the bandwidth of the controller has a magnitude of 10 dB. Therefore the 
controller will amplify input uncertainties at frequencies where flutter is most likely to occur. 

Robustness to additive uncertainty or gust responses is strongly affected by the peak of the closed-loop transfer 
function, KS 0 . Over the entire bandwidth, the singular values of KS 0 are below 0 dB and, thus, gust-like 
disturbances will be attenuated. The magnitude of the singular values of KS 0 fall off quickly after 6 rad/s. This 
predicts that higher-frequency turbulence will be rejected. 

Robustness to sensor uncertainty is predominantly dependent on the peaks of S 0 and T 0 . From 0.1 rad/s to 
1 rad/s the singular values of S 0 show that plant uncertainty is neither amplified nor reduced; however, near flutter 
frequencies the plant uncertainty is amplified. The singular values of T 0 indicate that sensor noise or output 
uncertainty will not be strongly amplified past 25 rad/s. 

Some aspects related to tracking may also be identified from Fig. 9. For best tracking accuracy, T 0 must be equal 
to 0 dB over most of the controller bandwidth. Figure 9 shows that there are several peaks above 0 dB, which may 
lead to overshoots or poor performance. The tracking history can be investigated by inputting doublet reference 
inputs to the controller. The set of charts presented in Figs. 10-12 expands on the performance and robustness 
characteristics of the controller. 
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Figure 10. Performance chart for doublet inputs to tracked variables: a) control surfaces; b) pitch angle 
tracking; c) bank angle tracking; d) SW1B modal displacement tracking; and e) SW1T modal displacement 
tracking. 
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Figure 11. Mu analysis over speed range: a) n chart; and b) corresponding uncertainty weightings. 
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Figure 12. Performance over 3.3% variation of velocity and uncertainty range: a) pitch angle tracking; b) 
bank angle tracking; c) SW1B modal displacement tracking; and d) SW1T modal displacement tracking. 
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It was desired that the responses of the system to a doublet have low rise time and small overshoots. This goal 
was for the most part achieved (see Figs. 10 and 12). All of the signals had low rise times of 0.25 s to 1.5 s. The 
change in pitch angle, 89, had excellent tracking characteristics. The change in bank angle, 8(f), showed a reasonable 
response. The tracked modal displacements were non-minimum phase. The modal displacement, 8SW1T 
experienced a 25 percent overshoot, which was predicted by the singular values of T 0 . 

There was very little rigid-body longitudinal and lateral coupling. It was interesting that the commanded 
rigid-body signals moved off their reference point when modes were commanded. The converse also happened; 
pitch angle commands tended to disturb the modal displacements. This indicates that a strong coupling between 
rigid-body motion and flexible motion in the models exists, which cannot be easily avoided. 

The normalized control movements for 89 and 8cp tracking were reasonably small. In fact, the movement of the 
control surfaces for 8(f) tracking was almost unnoticeable. During modal displacement commands, the normalized 
control movements moved to a maximum amplitude of 0.4 on the body flaps. Since the outputs are scaled to half of 
the actuator limits, this is a relatively large control movement, but is not unreasonable. In Fig. 10(a) it can be seen 
that the normalized control movement was significantly larger when commanding the modal displacements, than 
when commanding the flight variables. Significant control authority is required to perturb the X-56A structure from 
equilibrium. However, improper scaling may be to blame. Based on this analysis, the nominal performance (NP) 
appears to be adequate, although it might be improved with better weights and plant scaling. 

Figure 9 indicates that some disturbances and uncertainties are amplified; however, the actual impact of these 
uncertainties on the stability must be checked with a fi analysis as shown in Fig. 11. The fi analysis calculates 
robustness of the modal controller across dimensions of speed and uncertainty. For a range of plants G v additive, 
multiplicative input and output uncertainties with 10 percent weightings (see Fig. ll[b]) was added. The chosen 
shapes of the weightings were meant to replicate worst-case scenarios. Additive uncertainties may be more likely to 
occur in the 0.1 rad/s to 10 rad/s range due to gust disturbances. Multiplicative uncertainties may be more likely to 
occur at higher frequencies due to actuator dynamics and sensor uncertainties and noise. For the generalized plant 
with uncertainties, fi were then calculated across the operational frequency range with one controller, as shown in 
Fig. 11(b). 

Recall that RS is guaranteed if the condition given in Eq. (8) is met. Therefore, it is desirable that fi be less than 
one for all structured uncertainties. The maximum fi for all perturbed plants was 0.92. Therefore, RS is achieved. 
The variation of fi across the frequency range shows that it tended to increase at closed-loop flutter mode 
frequencies. The fi tended to increase for speeds past the design speed. Conversely, the \i decreased for plants below 
the design speed. This happens because as speed increases the damping of the flutter modes decreases (see Fig. 
4[b]). The reverse happens when decreasing in speed. Figure 11 indicates that the design requirement to be RS to a 
3 -percent variation in speed was met. Hence, it is shown here that the modal controller can be designed to be a 
robust to modeled uncertainties. Performance at off-design conditions must also be analyzed. 

Robust performance (RP) is achieved if and only if the performance objectives are achieved for all possible 
plants about the nominal plant up to the worst-case uncertainty. 19 There were no hard performance requirements 
defined other than that the tracked signals have low rise times and small oscillations. Figure 12 indicates that the 
responses from the perturbed plants was tight around the acceptable nominal plant responses shown in Fig. 10. The 
overshoots on the modal displacements for some perturbed plants did increase to 50 percent. But this overshoot 
increase is acceptable for the current study. Therefore, RP is achieved for these lenient performance requirements. 
Since NP, RP, NS, and RS are adequately achieved, the modal controller is acceptable. It meets basic criteria 
necessary for utilization in an aircraft. This completes Phase I. The controller designed here is used in the below 
discussion for virtual deformation control simulation. 

V. Virtual Deformation Control Simulation 

The modal controller derived above adequately satisfies the robustness and performance objectives for the 
X-56A model. The simulation study presented here demonstrates the use of the modal filter for virtual deformation 
control of the aircraft. To complete Phase II, the simulated fiber optic sensors are placed on the X-56A FEM. The 
sensor strain modes are calculated and the sensor strain mode matrix is formed. It is then shown that in order to 
create a proper deformation signal for Phase III, rigid-body translations and rotations must be removed from the 
deformation mode shapes. An objective function is given, which is minimized to achieve this goal. For 
demonstration, the controller designed in Phase I is put into a simulation with the modal filter designed in Phase II. 
The wing tips are commanded to a particular deformation using the reference signal derived in Phase III. 
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A. Fiber Optic Sensor Placement 

Traditionally, strain sensors or accelerometers used for active structural control are placed in an optimal sense. 
The FOS can measure strain at locations that are less than a one-half-inch apart, and can contain several thousand 
measurement stations along a single fiber. One of the benefits of having so much spatial sensor resolution is that 
optimal sensor placement (OSP) routines 32,30 lose significance. The first few structural modes in bending and torsion 
are most significant to flutter, so the sensors are placed from root to wing tip to capture bending information. The 
sensors are laid from trailing edge to leading edge to capture torsional motion. The resulting sensor configuration 
looks like three claw marks across each wing. The resulting sensor configuration is given in Fig. 13. 
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Figure 13. The fiber optic sensor layout on the X-56A model. 

With the six fibers placed, points on the aircraft were also chosen where deformation could be commanded with 
the reference signal derived in Phase III. These points can be anywhere on the aircraft (see Ref. 15). A multitude of 
points could have been chosen, but for simplicity four points were selected. While not required, an OSP technique, 
known as Kammer’s effective independence (El) procedure, 33 was utilized to locate the four points. This procedure 
starts by forming a sensor projection matrix from a selection of mode shapes. In this case, the sensor deformation 
mode shapes were utilized. The rows of the modal matrix were iteratively reduced based on the smallest diagonal of 
the projection matrix. After implementing this procedure, the resulting four deformation command locations were 
determined and are shown in Fig. 13. 

These points are right-wing trailing edge (RWTE), right-wing leading edge (RWLE), left-wing trailing edge 
(LWTE), and left-wing leading edge (LWLE). Their index locations in the modal matrix are referred to as l r . As in 
previous work 15 results indicate that the wing tips carry the most modal information. Note that the selected points do 
not have to be on sensor locations for the virtual deformation concept to be employed. 

B. X-56A Sensor Strain Modes 

The strain and deformation modes defined at the SFOS locations are necessary to develop an operational modal 
filter. To model sensor locations on the model, a spline was made on existing nodes on the FEM in MSC Nastran 
from root to tip on the wing. The distance between each node on the spline was set at one-half-inch, which was the 
spacing between each measurement grating on the fiber. 8 At each node of the spline or sensor location a minute 
mass was added. The sensor nodes on the spline were then interpolated with RBE3 MSC Nastran cards. The modal 
analysis conducted in MSC Nastran gives sensor deformation mode shapes corresponding to each flexible mode. 
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The close proximity of the sensors (of approximately one-half-inch) along the SFOS gives a unique opportunity 
to calculate axial strain modes directly from the modal deformation at the sensor locations. To our knowledge, this is 
not a common practice. It is also possible to calculate the strain-displacement relation using the finite element 
method 34 which may be more reliable for complex structures. One might also perform a static strain analysis for 
each deformation mode shape to obtain a strain mode shape. For simplicity, the following technique is utilized to 
convert each sensor deformation mode shape from MSC Nastran to a strain mode shape. Recall the axial strain in an 
element, 35 where x t is in the direction of the strain measure, shown in Eq. (14): 


du 

dx t 


(14) 


The FOS measure axial strain, so it is assumed that x t is in the local coordinate direction of the SFOS layout. Figure 
14 is presented as a visual aid for the breakdown of Eq. (14) into a usable format. 



Figure 14. The fiber optic sensor locations deformed for notional mode shape. 

Figure 14 represents sensors from the SFOS along the aircraft wing, which have been labeled for clarity. Each 
sensor has a line drawn between it, where a fictitious isotropic beam is assumed to exist. The deformed sensors are 
connected with blue lines. The original sensors are connected with black lines. Consider the deformed fictitious 
beam, to be in blue. Linear algebra may be used to show that the strain in any fictitious beam element of the SFOS 
shown in Fig. 14 is calculated as shown in Eq. (15): 


ProJutoPpP'-QP 

OP 


where OF is a directed line segment from the un-deformed sensor location OeM? to the un-deformed sensor location 

PeM 3 , and 0 r P r is a directed line segment from the deformed sensor location O'eM 3 to the deformed sensor location 
P'eM 3 and u is the unit vector. The projection in Fig. 14 calculates the final length of the beam section with respect 
to its original beam orientation. The difference between the original length and the deformed length, divided by the 
original beam length, gives the strain in the beam element. The strain in one beam may differ from the strain in the 
beam next to it. The strain from modal deformation on the r th sensor (or sensor location 0 ) may be defined to be the 
average of the strain on either side (see Fig. 14) as shown in Eq. (16): 
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where e xx is the axial strain measured on the fore side of the sensor and £ xx is the axial strain measured on the aft 
side of the sensor. This operation can be carried out for each mode and at each sensor node, and strain modes x/j m 
may be formed by repeating Eq. (16) for each sensor. The sensor strain modes that result are collected into the 
desired strain mode matrix defined at SFOS sensor locations, as shown in Eq. (17): 

VfOS (17) 


To compute the strain at point 0 with good precision, candidate sensor locations P+ and Pi must be selected close to 
the original sensor. Since the change in strain was so small between two adjacent sensors, two sensors were selected 
a few inches away from the desired sensor location. Even with this technique, the resulting strain modes from Eq. 
(16) appeared noisy. Thus, the strain was fit to a polynomial surface for each wing and an M-estimator 36 was utilized 
to trim outliers. The first four sensor strain modes used for the X-56A simulation studies are given in Fig. 15. 
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Figure 15. Sensor strain mode shapes: a) SW1B; b) AW1B; c) SW1T; and d) AW1T. 

The sensor strain modes were computed for all 14 flexible modes of the X-56A, however, only the first four 
modes are shown here. The SW1B mode has a strain distribution which starts out very high near the root of the wing 
(see Fig. 15) and goes to nearly zero at the wing tips. It is typical that in bending, strain is highest near the root of 
the beam. 16 The AW IB mode shows an antisymmetric strain distribution in bending. The SW1T mode shows a 
distribution of strain which has an inflection near zero close to the mid-span of the wing, indicating that for torsion 
modes where the wing is twisting, very little strain exists in the mid-span of the wing. The AW IT mode, like the 
AW IB, mode has an antisymmetric strain distribution in torsion. The sensor modes are collected into the modal 
matrix defined in Eq. (17) and, together with Eq. (6), the modal states of the aircraft can be measured. With the 
strain mode matrix prepared, attention is turned to the reference signal, which is needed for virtual deformation 
control. For reference signal creation, the deformation mode shapes shown in Fig. 3 must be transformed. 

C. Reference Signal Creation 

Phase III from the modal filter control design methodology is to design a modal reference command. From Eq. 
(5) a least-squares approximation may be used and the modal reference command is shown in Eq. (18) (see Ref. 15): 
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Qref = (0 7 (/ r ,m r )0(/ r ,m r )) 1 <P T (l r ,m r )d re f 


( 18 ) 


where d re f is a vector of deformations from the un-deformed aircraft corresponding to the index vector, l r , in the 
pure elastic deformation modal matrix, O. Recall that each row of the modal matrix corresponds to a physical 
location on the aircraft, (x c ,y c ,z c ). The index of modes, m r , corresponds to the index of modes within the modal 
matrix, which the controller is designed to track. In this case, m r corresponds to {1,3}, since the SW1B and SW1T 
modal displacements are tracked. 

Performing virtual deformation control on the X-56A model is not straightforward because the flexible modes of 
the vehicle are not purely elastic. The normal mode shapes (see Fig. 3) are corrupted by strong rigid-body motion. 
The use of these modes to recreate the modal reference signal will lead to large displacement errors. The mode 
shapes required in Eq. (18) are calculated by an optimization procedure on the original mode shapes. The 
optimization procedure rotates and translates every node in the mode shape about the aircraft center of gravity (eg). 

Recall that a mode shape corresponds to deformations at N grid points. At each grid point, the mode shape is 
defined for six DOF, so each mode shape technically has 6 N DOF. The number of DOF are increased with global 
DOF parameters so the total DOF in each mode equate to 6N + 6 DOF. The global DOF are the only variables in 
the objective function defined for the i th mode shape given by Eq. (19): 
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where x cg , y cg and z cg are the scalar coordinates of the eg of the aircraft, T(-) is a right-hand-rotation matrix 
operator, x U9 y u and z u are N x 1 coordinate vectors of the un-deformed aircraft, x^, y^,, and are N x 1 
coordinate vectors of the i th modal deformation vector, The six scalar variables in the optimization function 
include translational coordinates: 8x c , 8y c , and 8z c , and rotations about all three axes: 86 c , 8 <p c , and 8i/j c . 

By minimizing the objective function (see Eq. [19]) for the i th mode, with respect to these six scalar variables the 
relative distance between the un-deformed aircraft and the deformed aircraft is minimized. Since the squared 
Euclidean norm is a strictly convex function of its respective arguments, the objective function is also convex; 37 
therefore, there are numerous discrete optimization methods available in the literature which can be used to solve it. 
For each mode i, the objective function is minimized with a locally convergent random search technique known as 
the compass method. 37 Four of the resulting corrected mode shapes are given in Fig. 16. 
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Figure 16. Corrected mode shapes a) SW1B; b) AW1B; c) SW1T; and d) A WIT. 

The mode shapes with rigid-body components removed (see Fig. 16) may be used with Eq. (18) to develop 
appropriate reference signals. Note that the flexible movement within each mode appears amplified when compared 
to Fig. 3. Deformation commands can now be adequately transformed to achieve virtual deformation control. 

D. Simulation 

The simulation architecture introduced in Fig. 6 represents the inner-loop design for the X-56A simulation model. 
The inputs to the controller are references, provided by a notional outer-loop control system. The references in Fig. 6 
are separated into rigid-body commands, x r A J , from the flight computer and vertical deformation commands, d re f, 
at selected locations (see Fig. 13), namely, LWLE, LWTE, RWLE, and RWTE. The reference vertical deformations 
were chosen to be 0.6 percent of the span of the vehicle, which corresponds to modal displacement commands of 
approximately twice the scaled value. This value was thought to be reasonably small, since the control movement 
for scaled modal displacement commands was less than 50 percent of the actuator range as shown in Fig 10(a). 

Since the reference deformations were equal, a bending-type motion for the aircraft wing shape results. Different 
deformation choices could lead to torsion-type motion instead. The deformation references can be verified to 
generate a bending motion when transformed with Eq. (18) into a modal reference. For the following study, the 
airframe state references are set to 0. 

The controller that was designed above (see Figs. 9-12) takes the difference of the references and the feedback 
signals to produce a control input signal. Multiplicative noise is added to the controller inputs, with a mean of 0 and 
standard deviation of 0.1. The plant takes the inputs via Eq. 3 and produces the state vector, x q , in the output, which 
includes the measurements SW1B, SW1T, and the residual modes. The plant also outputs the rigid-body 
measurement vector, x AF . 

The SFOS measurements are simulated by multiplying the sensor strain matrix, , with the vector, x e . 

Normal noise, n s , with mean 0 and a standard deviation of 3 fi s is then added to this strain vector. Since the rigid 
measurements are scaled, a normal distribution with mean 0 and standard deviation of .01 is added to the 
measurements to simulate noise. The results of the simulation to a reference doublet command are shown in Fig. 17. 
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Figure 17. Virtual deformation tracking time history: a) deformation tracking at wing tips; b) tracked modal 
displacements; and c) airframe states. 
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The change in deformation near the wing tips (see Fig. 17 [a]) indicates that the actual deformations moved very 
close to the deformation set-points. The error is due to the modal displacements from residual modes (see Ref. 15). 
This is expected, since the trailing-edge actuators are not collocated with the structure’s sensors. Due to this fact and 
since the modes are coupled aerodynamically, independent modal control cannot be achieved. 16 

The commanded modal displacements are shown in Fig. 17(b), and they are tracked well. Since the command is 
primarily due to bending, the SW1T modal reference is very small. The scaled SW1B modal displacement 
amplitude moves the greatest amount. Both modes had some overshoot and oscillations during the initial part of the 
command, which was predicted and is seen in in Fig. 10. The steady state response has oscillations primarily due to 
the rather large multiplicative input noise modeled in the system. 

During the shape-deformation maneuver, the aircraft experienced no enduring change of velocity or pitch angle. 
The angle of attack, however, achieved a steady-state change; the change in angle of attack should have also 
effected an enduring change in velocity. An angle-of-attack change generally leads to a drag increase or decrease, 
however, this was not observed, likely due to improper drag modeling. If the results were taken as truth, it could be 
concluded from this analytical model that the effective lift over drag (L/D) ratio has changed. This observation could 
be a result of control surface movements, the structural deformation itself, or both. The caveat here is that the results 
are achieved on a linearized plant model defined for small perturbations around the trim point. In flight, the aircraft 
achieves a new trim point from the angle-of-attack change, which may result in a lower or higher L/D ratio. 

Outer-loop control requirements are also observed from Fig. 17. The results indicate that the flight-path angle 
has changed due to shape deformation commands. Therefore, the aircraft will require a heading-angle-hold 
outer-loop controller to make virtual deformation control feasible. The outer-loop controller could command the 
pitch angle to control the flight-path angle; or it could command engine throttle settings, which are not currently 
modeled. If the throttle is changed to achieve trim, fuel efficiency degradation or improvement may be realizable, 
depending on the command to the structure. 

Previously, it was shown that the X-56A model experiences strong flutter at the design speed (see Fig. 5). 
Technically, the shape controller demonstrated above is capable of suppressing flutter. From the previous plots, this 
has not been shown explicitly. To demonstrate this in a traditional sense, the aircraft is perturbed as before and the 
control system is turned on after a period of time of 1.6 s. To reduce control-induced oscillations, it was ramped in 
over a 1-s period, with all references set to 0. In 2 s, the shape tracking controller achieves AFS, as shown in Fig 18. 


28 

American Institute of Aeronautics and Astronautics 



scaled amplitude scaled amplitude scaled amplitude 


30 


5SW1B 



Figure 18 Active flutter suppression time history: a) BFF suppression; b) SWBT suppression; and c) AWBT 
suppression. 
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During AFS there was an altitude change (see Fig. 18 [a]); however, this change occurred because an outer-loop 
guidance system is not currently modeled. The achieved flutter suppression was not unexpected, since control was 
implemented on a linear plant model. Nonlinear effects have not been included. Only the poles must be considered 
in such a model. The poles of the closed-loop system have real negative parts and the closed-loop system is, thus, 
stable. The open-loop poles and the closed-loop poles are plotted together in Fig. 19. 




real part real part 

Figure 19. Open-loop versus closed-loop poles: a) open-loop poles; and b) closed-loop poles. 

The poles of the open-loop system indicate the presence of the three flutter modes identified above (see Fig. 4). 
There are also poles from lightly-damped structural modes near the vertical axis starting at 54 rad/s. When the loop 
is closed with the controller, the flutter poles migrate to the stable region in the left half-plane. The lightly-damping 
structural mode poles remain unchanged, as they are outside the controller bandwidth. The controller introduces 
poles seen tightly grouped near the lower right-hand corner of Fig. 19(b). Integrator and estimator poles are present. 
The lowest damping of the closed-loop system poles within the actuator bandwidth is 0.08. This represents a 
significant improvement over the open loop BFF damping of -0.13. 

The results of all simulations must be accepted as notional. Recall that all of the simulations are completed on a 
linear plant model with simulated strain measurements. Further work with time-marching nonlinear simulations will 
be necessary to validate the controller architecture for AFS and active shape control, however, these preliminary 
simulation results are promising. Time histories indicate that the modal filter is a suitable estimator for both AFS 
and virtual deformation control. Simulation results similar to this were published previously in Ref. 15, but those 
results were for a clamped-wing model only. 

VI. Conclusions and Recommendations 

There were three primary objectives set out to be accomplished in this paper. The first objective was to 
demonstrate that the least-squares modal filter is a candidate aircraft modal estimator and that it can be adequately 
used for control feedback. For practical implementation, three phases for incorporating the modal filter were 
presented and followed throughout the paper. The modal filter was demonstrated with simulated fiber optic sensors. 
The estimated modal displacements were used for both shape control and active flutter suppression of the X-56A 
simulation model. It was found that significant angle-of-attack changes occurred during shape control. The angle-of- 
attack change may lead to improved lift to drag ratios on the aircraft. Therefore, structural shape control may 
someday enhance performance in fuel efficiency or maneuvers. The caveat is that it is not known what lift to drag 
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ratio will result at the new trim condition. Nonlinear simulations and experimental studies are required to explore 
this further. 

The second objective was to demonstrate that the modal filter and distributed sensing can be used to control the 
deformation of the aircraft wing structure in flight. By designing a controller for the modal filter and simulated fiber 
optic sensors, and commanding the wing-tip deflections, this objective was satisfied. Only bending commands were 
presented here, however. The deteriorating effects of residual modes may be more prominent with other command 
shapes. Modeling errors could also appear prominently in a realistic bench case study. In-flight experiments are also 
required to validate its use for aircraft state estimation and control. 

The third objective was to show that the modal controller can be designed to have good performance and be 
robust to parameter variations. This objective was satisfied by employing the powerful mu synthesis technique. It 
was shown that the controller had robust performance and was robustly stable to structured uncertainty and model 
speed variations. The requirements for robust performance were lenient, however, and further improvements can 
certainly be made here. Future endeavors will extend the controller to be robust to fiber optic sensor failures. 
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